Phase separation in an homogeneous shear flow: 
Morphology, growth laws and dynamic scaling 



Ludovic Berthier 
Laboratoire de Physique, Ecole Normale Superieure, 
46 Allee d'ltalie, F-69007 Lyon France and 
Departement de Physique des Materiaux, U. C. B. Lyon 1, 
43 Bd du 11 Novembre 1918, 69622 Villeurbanne Cedex, France 
(February 1, 2008) 

We investigate numerically the influence of an homogeneous shear flow on the spinodal decompo- 
sition of a binary mixture by solving the Cahn-Hilliard equation in a two-dimensional geometry. 
Several aspects of this much studied problem are clarified. Our numerical data show unambiguously 
that, in the shear flow, the domains have on average an elliptic shape. The time evolution of the 
three parameters describing this ellipse are obtained for a wide range of shear rates. For the low- 
est shear rates investigated, we find the growth laws for the two principal axis R±{t) ~ constant, 
7i|l (t) ~ t, while the mean orientation of the domains with respect to the fiow is inversely propor- 
tional to the strain. This implies that when hydrodynamics is neglected a shear flow does not stop 
the domain growth process. We investigate also the possibility of dynamic scaling, and show that 
only a non trivial form of scaling holds, as predicted by a recent analytical approach to the case of 
a non-conserved order parameter. We show that a simple physical argument may account for these 
results. 
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I. INTRODUCTION 

The study of phase ordering kinetics has a long his- 
tory . The canonical example is the coarsening pro- 
cess follovifing the quench of a binary mixture A-B be- 
low its spinodal line. The properties of the resulting do- 
main growth are rather well understood . In the case 
where A and B are in equal concentration, an isotropic 
bicontinuous structure emerges, which is characterized by 
a typical length scale L(t) growing as a power law of time 
t. Moreover, as t increases, the structure evolves in a self- 
similar manner in the sense that its statistical properties 
are the same when space is rescaled by L{t). Any function 
C{r,t) that depends on space and time is then a function 
the reduced variable \r\/L{t) only: C{r,t) = C(\r\/L{t)). 
This property is termed 'dynamic scaling' |l|-|2]- 

The study of the spinodal decomposition in an ho- 
mogeneous shear flow is of fundamental and practical 
interest 0], but despite an enormous amount of exper- 
imental KI numerical |p^-p3|| and analytical works 
|g,|2|j2|-g^, the problem is still not fully settled_p6|[ 
Early experiments pJ§;0| and simulations ||lT 



have shown that isotropy is lost, the morphology of the 
bicontinuous structure being elongated along the flow di- 
rection. Hence, a single length scale can not describe the 
full structure. There are thus several questions which 
naturally arise. (1) How to characterize quantitatively 
the growing structure? (2) What is the time dependence 
of the different length scales, as compared to the un- 
sheared case? (3) Does the growth stop after some time? 
(4) Does a suitable generalization of the dynamic scaling 
property hold? 



It is interesting to remark that in spite of the large 
number of works cited above, a definite answer to all 
these questions is still lacking. Several reasons make this 
problem non-trivial. First, the role of hydrodynamics is 
far from being understood. Simple physical arguments Q| 
show that it may become preponderant at large times, 
predicting a saturation of the domains at a typical size. 
This assumption is resolved neither at the experimental 
level (see, e.g., the opposite conclusions of Ref. and 
Refs. nor at the numerical level, where powerful 

algorithms are needed to correctly deal with hydrody- 
namics. At present, this limits the numerical analysis to 
too small sizes to make any definite answers to questions 
(l)-(4) above, although much progress has been made 
recently 

Second, at the theoretical level, no analytical solution 
of a reasonable model of spinodal decomposition (even 
neglecting hydrodynamics) is available. One has then to 
make some predictions from the solution of solvable, but 
less realistic models, like the 0{n) model in the large-rt 
limit p^ , or from the approximate solution in the case 
of a non-conserved order parameter [^5| . A scaling argu- 
ment based on the hypothesis that a generalization of dy- 
namic scaling holds has also been developed in Ref. [|2| . 

Third, there are, to our knowledge, no numerical sim- 
ulations (neglecting hydrodynamics) validating these an- 
alytical predictions. Moreover, the scaling hypothesis on 
which is based the analytical argument of Ref. ||l^ has 
only been tested in Ref. ||], with negative results. Thus, 
the validity of the predicted growth laws may also be 
questioned. More crucially, there is up to now no con- 
sensus concerning the numerically measured growth laws: 
we discuss this point in a more detailed way in the last 
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Section. 

In this work, we numerically study the spinodal de- 
composition process in a shear flow by solving the Cahn- 
Hilliard equation in two dimensions. All hydrodynamic 
effects are neglected. Although this involves a drastic re- 
duction of the experimental situation, it is necessary, in 
our opinion, to have a good understanding of this 'sim- 
ple' case, before studying more realistic problems. Our 
algorithm is different of but our study is technically com- 
parable to the most recent one We shall however 
explore a wider range of shear rates, and this will lead us 
to a different interpretation of the numerical data. 

The paper is organized as follows. In the next Section, 
we define the model and describe the numerical proce- 
dure to solve it. Section III briefly recalls the results 
obtained when the shear flow is absent. Section IV de- 
scribes the morphology of the domains under shear and 
its time evolution. Section V focuses on the problem of 
dynamic scaling. In the last Section, we compare our re- 
sults with the relevant existing data in the literature and 
give a simple physical argument to explain them. 



II. MODEL AND DETAILS OF THE 
SIMULATION 

In this work we focus on the standard model for the 
spinodal decomposition of binary mixtures and solve nu- 
merically the Cahn-Hilliard equation 



dt 



5F[<f>] 
50(r,t) 



(1) 



In this expression, the order parameter 4>{r, t) is a scalar 
quantity which can be linked to the concentration ca {cb) 
of the component A (B) of the mixture by the relation 
(j) = \ — 2c A = "^cb ^ ^- Eq. (|l|) has the form of a con- 
tinuity equation, which implies that the order parameter 
is a conserved quantity. The free energy is of the 
Ginzburg-Landau type 
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where the equilibrium correlation length ^ has been in- 
troduced. The noise term ry is a random Gaussian vari- 
able, characterized the two moments (77 (r,t)) = and 
{r]{r,t)r]{r',t')) = -2TS{t~t')V^5{r-r'); T is the tem- 
perature of the thermal bath. All the simulations will 
be carried out at T = 0, since temperature is essentially 
irrelevant in this process ||^ . It is also known to delay the 
onset of the asymptotic regime . The second term in 
the left hand side of Eq. (|l|) results from the advection of 
the order parameter by the velocity field. The case of an 
homogeneously sheared system will be investigated. The 
flow is taken to be in the a;-direction and the velocity 
field is then v = jy Bx, which defines the shear rate 7. 
We shall concentrate on the case of a constant shear rate. 



We want then to solve numerically in two spatial di- 
mensions the following equation: 



dt 



-ly 



dx 



- rv^ (c^v^ 



(3) 



where both space and time dependences have been re- 
moved for clarity. This is done by combining the nu- 
merical methods of Refs. A new frame [x' ,y') 
is first defined by |^ x' = x — S{t)y, y' = y, where 
S{t) = di'7 is the strain. In the case of a constant 
shear rate, it is simply given by S{t) = ^t. Defining 
further (t){r,t) = ${r',t), Eq. (|) becomes |2| 



dt 



-V^ V 



(4) 



with 



_a_ _9_ _d_\ 



_d_ 



— -S(t) — 
dy' ^'dx' 



(5) 



After transformation, Eq. (Q) is formally identical to 
the Cahn-Hilliard equation without shear, which is solved 
by the implicit spectral algorithm developed in Rcf. . 
Space is measured in units of the correlation length ^ 
(also the interface width) and time in units of ^^/F. This 
microscopic time scale represents the typical time it takes 
to create a well-defined domain wall. Periodic boundary 
conditions are imposed in the deformed frame. The sin- 
gle parameter of the simulation is then the shear rate 
7, which introduces a time scale 7"^. The choice of the 
parameters for the discretization has been discussed in 
Rcf. (2|] and the values At = 0.5, Ax = Ay = 0.5 are 
used throughout the simulation. For each shear rate, the 
system size has been carefully checked to be large enough 
so that the reported growth laws are unaffected by the 
boundaries. Since the growth is strongly anisotropic, a 
rectangular simulation box has been chosen with sizes up 
to Ly = 512 and = 8192. The shear rates investigated 
in this paper are 7 = 0.04, 0.02, 0.01, 0.005, 0.0025 and 
0.00125. This corresponds to a time scale 7^^ in the 
range [25, 800]. We wish to emphasize that the condition 
7~^ ^ C^/r has to be fulfilled, since we are interested 
in the scaling regime where well-defined domains coarsen. 
This remark will become important for the interpretation 
of the numerical results. An alternative solution would 
be to apply the shear flow after an initial transient, so 
that large domains have already grown. Such initial con- 
ditions are discussed at the end of the paper. 



III. ZERO SHEAR CASE 

Although the 7 = case has been extensively studied 
||l|-^, we briefly consider this well-known situation with 
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three aims. These results are presented to (i) vahdate our 
numerical procedure, (ii) present the quantities of inter- 
est and above all (iii) make comparisons to the sheared 
case easier. 




FIG. 1. Snapshots of size 256 x 256 of the unsheared spin- 
odal decomposition at times t — 70, 381 and 3761 (from left 
to right). Each color represents one phase of the mixture. 

The domain growth takes place typically as in Fig. ||, 
where an isotropic bicontinuous structure coarsens with 
time. This coarsening process is well characterized by 
the two-point correlation function defined by 

1 

V 



C{r,t) 



d'^x{(f>ix,t)(f){x + r,t)), 



(6) 



which is nothing but the Fourier transform of the struc- 
ture factor, experimentally measured through light scat- 
tering experiments. A typical two-point function is repre- 
sented in Fig. 1^ which shows the isotropy of this surface. 
The average shape of the domains of Fig. |^ may be ex- 
tracted from this plot by taking the intersection of this 
surface with a horizontal plane z = constant. This al- 
lows to measure the time dependence of the length scale 
Rx{t) {Ry{t)) in direction x (y). Both length scales are 
represented in the inset of Fig. |[ and have the expected 
power law behavior ~ Ry oc t-^^^ [0^- 

The dynamic scaling hypothesis is tested in the main 
frame of Fig ^, where the two-point function is circu- 
larly averaged and plotted as a function of |r|/L(t). This 
works perfectly well. 

We are thus confident in our numerical setup, and we 
shall now address the question of the influence of the 
shear flow on the spinodal decomposition. 





FIG. 2. Two-points correlation function, Eq. (^), in the 
7 = case, at time t = 100. 



FIG. 3. MAIN: Circularly averaged two-point correlation 
as a function of the rescaled variable \r\/L{t) for differ- 
ent times in the range [50,1000]. By construction, one has 
C{L{t),t) = 0.2, the choice of 0.2 is indifferent. INSET: 
Growth laws in the x (circles) and y (triangles) directions. 
The dashed line is a fit to a power law t^^^ . 



IV. ANISOTROPIC GROWTH OF THE 
DOMAINS UNDER SHEAR 

A. Basic observations 

The time evolution of the domains when a shear flow 
is applied is followed in Fig ^. This evolution is the basic 
result of all previous numerical works p^-|2^] and shows 
that the domain growth is essentially unaffected for times 
t < since the first snapshot is very similar to those in 
Fig. |l|. At intermediate times, the domains begin to have 
an anisotropic shape: an average direction of the domains 
is clearly apparent. This direction rotates and becomes 
more aligned with the flow when the strain increases. At 
large times, t ^ 7"""^, domains are nearly aligned with 
the flow, and have a strongly anisotropic shape. 

These features are also clearly discernible in Fig. ||, 
where we have plotted the two-point correlation function 
for a strain S{t) = 5 at 7 = 0.01. The surface is clear! 
stretched in the flow direction, compare with the Fig 
In particular, it becomes impossible to perform a circular 
average as in the unsheared case: isotropy is lost. 

The average shape of the domains is again recorded 
through the intersection of the two-point correlation 
function with a horizontal plane. For definiteness, we 
take z = 0.5, but the value 0.5 is unimportant. It is 
found numerically that this intersection is very well rep- 
resented by an ellipse. Computing the parameters of this 
ellipse gives then access to two typical length scales, (t) 
(large axis) and R± (t) (small axis) , and to the mean ori- 
entation of the domains, 9{t) (angle between the large 
axis and the a;-direction) . We present the typical time 
evolution of the elliptic shape of the domains in Fig. 0, 
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where data are also fitted to an elliptic form, in a very 
satisfactory way. 

These three quantities depend both on the time t, and 
on the shear rate 7 and we successively study in the fol- 
lowing subsections the time evolution of 9{t) and of the 
length scales R\\{t) and R±{t). 




FIG. 4. Snapshots of size Ly = 512 and Lx = 2048 (parts 
of a 512 X 4096 system) for a shear rate 7 = 0.01 and strains 
S(t) = 1, 5, 10 and 50 (from top to bottom). Each color 
represents one phase of the mixture. 



C{r,t) 




FIG. 5. Two-point correlation function for 7 = 0.01 and 
S{t) = 5. The surface is stretched in the a;-direction. Note, in 
particular, that the x and y ranges are different in this figure. 



Before performing this analysis, a remark has to be 
made about the identification of the relevant length 
scales. It is quite clear from the above analysis that 3 



parameters are needed to fully characterize the growing 
structure under shear. This feature was already noted 
in experiments 1^,^. In some of the previous numer- 
ical works, only two parameters were studied, namely 
the length scales in the x and y directions. Although 
these length scales should well represent the structure at 
long times, it is physically preferable to study instead 
the domain size in the directions defined by the angle 
9{t). In Ref. the direction of the domains' shape 
was recorded, but there was no attempt at a quantitative 
analysis of its behavior. We note last that this quantita- 
tive analysis of the domain morphology in two dimensions 
naturally arises from the analytical work of Ref. , in 
the case of a non-conserved order parameter. 
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FIG. 6. Time evolution of the intersection of the two-point 
correlation with a horizontal plane, for strains S{t) — 5 
(squares), 10 (circles) and 20 (triangles) and shear rate 
7 = 0.01. The points are the data, while the lines are fits 
to an elliptic shape. Note, in particular, that the x and y 
ranges are different in this figure. 



B. The mean orientation 0{t) 

The first effect of the shear flow is to give the domains 
an anisotropic shape and hence to create a preferred di- 
rection in the system. Physically, it can be expected 
that for a given shear rate 7, the mean orientation of the 
domains is inversely proportional to the strain S{t), as 
would be the case for a rigid rod advected by the shear 
flow: 

Oit) - (7) 
7t 

This relation is tested in Fig. |^, where ^0(7) is used as 
an adjustable parameter. This figure shows that the rela- 
tion (|^) is well satisfied in the whole range of shear rates 
investigated. 

The parameter is found to be a slowly increasing 
function of the shear rate. One finds 0o (0.00125) ~ 0.975 
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and 00 (0.04) ~ 2.0. By definition, this angle corresponds 
to the mean orientation of the domains when the strain is 
1, 9o = 6{t — 7~^), and its variation may be understood 
by the following argument. For times t < 7~^, the do- 
main growth has been mainly unaffected by the flow and 
thus the larger the shear rate, the smaller the domains 
at time t 7~^. Since large domains are more easily de- 
formed than small ones [because of the surface tension] , it 
is expected that at strain S{t) ~ 1, large domains (small 
7) are more deformed than small ones (large 7). Hence, 
the smaller the shear rate, the smaller 9q. 



a± = has been compared to the data. Second, fol- 
lowing the available analytical results |12 2jj2^, we also 
tried to fit the data with the ansatz 
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FIG. 7. Mean orientation of the domains as a function of 
the strain S{t) = 'yt. The straight line is l/'yt. 



C. Growth laws 

We turn now to the time evolution of the two length 
scales i? j_ {t) and R\\{t). These quantities are studied for a 
very broad range of shear rates from 7 = 0.04 (7"^ = 25) 
to 7 = 0.00125 (7^^ = 800). All our results are summa- 
rized in Fig. H, where and are represented for 
each value of the shear rate as functions of the strain. 
We obtain numerically that the growth laws are well rep- 
resented at large strains by the following algebraic forms, 



[t) ~ i?||o ■ (7t)"" , R±{t) ^ R±o ■ iity 



(8) 



which define the growth exponents and a±. This 
was first obtained in a simulation by Padilla and 
Toxvaerd |22| and subsequently in similar systems in 
Refs. [|jl|l|,|2|. 

In FigT^, we have fitted our data using the forms (H). 
These fits deserve some comments. The growth law for 
becomes a nice straight line in a log-log plot at large 
strains S{t) > 10. Only the data points for such large 
strains have been used to compute the exponent ay . Con- 
cerning the time behavior of R± (t) , two different fits have 
been tested. First, a horizontal line corresponding to 



= -1. 



(9) 



Note that both fits are equivalent when a|| ~ 1, which is 
nearly the case for the two smallest shear rates. 
Several comments are in order. 

1 - The algebraic fits are clearly a very good repre- 
sentation of the data. Some logarithmic corrections that 
come out from the analytical study of the 0{n) model 
and of the non-conserved case may be present, but we do 
not expect to be able to determine them numerically. 

2 - Although the relation (^) reasonably accounts for 
the data, the value a± = is also possible and works 
well for all the shear rates investigated. This means that 
it could be possible to rescale all the curves for R±{t) 
by plotting R±{t)/R±o as a function of the strain for 
different shear rates, using R±q as a fitting parameter. 
This parameter is found to be a decreasing function of 
the shear rate, which means that the smaller the shear 
rate, the wider the domains. This rescaling is performed 
in Fig. ^, and the data are indeed compatible with this 
hypothesis. Let us add the remark that R±{t) varies (at 
most) by a factor 2 in all the simulations. This strongly 
supports the hypothesis that there is in fact no growth 
in the perpendicular direction. Of course, once again, 
logarithmic corrections can not be numerically dismissed. 

3 - We do not see any evidence for the oscillations re- 
ported in Refs. We do not have a clear explana- 
tion for that, but an hypothesis is that these oscillations 
are a preasymptotic artifact of the measure process it- 
self. The definition of the length scales used in Ref. jlj] 
leads in some cases to a ratio Rx{t)/Ry{t) ~ 10 for a 
strain S{t) ~ 1 For this strain, the domains are still 
nearly circular and in our simulations, this ratio is never 
above the value 1.5. However, this argument is made 
weaker by the observation that Qiu et al. |9| have used 
the same measurement procedure as Corbcri et al. [ p2[ , 
at a lower shear rate 7 — 0.01 and do not observe any 
oscillations. 

4 - An important point is the fact that the exponents 
apparently depend on the value of the shear rate. More 
precisely, we find that a\\ decreases from the value a\\ ~ 
1.35 at 7 = 0.04 to a\\ ~ 1.0 at 7 = 0.00125. This means 
that we are in fact measuring effective exponents and that 
the true asymptotic behavior has not been reached in 
some of the cases studied here. 

The problem is then to determine which exponent is 
the right one! As emphasized in the Introduction, well- 
defined interfaces exist only if 7 ^ 1 . This indicates that 
the 'true' asymptotic behavior is reached for the smallest 
shear rates investigated, and favors the value a\\ — \ 
found for 7 = 0.0025 and 0.00125. Together with the 
behavior of and the observation that the domains are 
wider at lower shear rates, we are led to the conclusion 
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that the growth exponents are given by 
a|| — 1, aj^ = 0. 



(10) 



The analysis of dynamic scahng properties in the follow- 
ing section will reinforce this conclusion. 
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FIG. 8. The two length scales R±{t) and R\\{t) are repre- 
sented as a function of the strain. Each figure is labeled by the 
corresponding shear rate. The symbols are the data, and the 
full lines are algebraic fits with the exponents indicated near 
each fit. For the strongest shear rates, a horizontal dashed 
line has been added as a fit to a± = 0. 
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FIG. 9. Test of the relation (|8|) for the smaU axis of the 
eUipse. The full line represents the case a± = 0, whereas the 
dashed line is for ax = 1/3. The symbols are the same as in 
Fig. 1^. The latter is clearly inconsistent with the numerical 
measurements for all the shear rates. 



V. DYNAMIC SCALING 

The last question we wish to address is the problem 
of dynamic scaling. We have recalled in Fig. (where 
7 = 0) that the two-point correlation function has the 
property that it can be rescaled in the form of a single 
variable function 



C{r,t)=C 



Lit) 



(11) 



This indicates that L{t) in the only relevant length scale 
in the asymptotic regime characterized by ^ <C L{t). This 
behavior is the basis of a scaling argument which, in 
the unsheared case, allows an elegant derivation of the 
growth laws |^ . 

In the shear flow, there are two relevant length scales 
and the scahng (11) can thus no longer be true. Two 
different generalizations have been suggested by analyt- 
ical works. The solution of the 0{n) model predicts the 
scaling form p3] 



Cir,t)=C 



y 



R,{ty Ry{t) 



(12) 



where Rx and Ry are the typical sizes in the directions x 
and y respectively. A different scaling is expected from 
the non-conserved order parameter case, namely ||25| 



C{r,t)=C 



R\\{ty R^{t) 



(13) 



where the subscripts refer to the rotating frame described 
previously. Corberi et al. jl^ have used the form (|l^) as 
a starting point to generalize the argument of Bray |3G] 



to the sheared case. It is thus important to see if this 
scaling behavior is detected in the simulation. Let us 
note that the tilt angle 9{t) is very small in the long 
time regime we are interested in. Then, one might ask if 
the difference between the two forms (p^-(p^ proposed 
above is relevant. Since the domains are very elongated in 
the x-direction, then even with a small angle there might 
be differences between the 'parallel' and the x-directions. 
Concerning the y and the 'perpendicular' directions, the 
support of the correlation function in these directions is 
very small, so that differences between the two are indeed 
not observable. 



1 

0.8 
0.6 







0.2 



-0.2 



' r- 


1 1 1 


1 1 1 






7( = 20 □ 


\ 


7 = 0.04 


7( = 40 


Q 


7* = 60 V - 






7* = 80 i 














1 


1 1 1 


1 1 1 



4 




xlR,{t) 

FIG. 10. Test of dynamic scaling in the a;-direction for 
7 = 0.04 and 0.0025. In both cases, the collapse is not satis- 
fying. 



We now present our numerical results. In Fig. |l^, we 
show for two different shear rates the attempt at the 
rescaling of the two-point function in the x-direction. In 
each case, we have considered a time window where the 
hypothetic scaling might hold {i.e. the growth law has 
reached its algebraic asymptotic form, S{t) > 10). We 
have also chosen Rx{t) in order to get the best collapse 
of the data. Clearly, the scaling ( [I4 ) does not work, in 
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both cases of a high and a low shear rate. 

We investigate the possibihty for the second scahng 
form to hold in the direction defined by the tilt angle 
0{t) in Fig. There is here a clear qualitative difference 
between the shear rate 7 — 0.0025 and 7 — 0.04, since 
the collapse of the data is excellent for the small shear 
rate, whereas there is a clear systematic evolution for the 
highest shear rate. 

In our opinion, these results are a very good indication 
that the scaling form (|2|) does not describe the asymp- 
totic behavior of the Cahn-Hilliard equation. Moreover, 
they show that a self-similar asymptotic regime has been 
reached for the lowest shear rates investigated, charac- 
terized by the scaling (|l^) of the two-point function. 

This is also confirmed by the inspection of the scaling 
properties in the 'parallel' direction (which is equivalent 
to the y-direction). In Fig. we show the results for 
this direction. Once again, the collapse is very good for 
the low shear rate 7 = 0.0025, whereas it is clearly not 
satisfying for a higher shear rate 7 — 0.04. 
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FIG. 12. Test of dynamic scaling in the parallel (or equiv- 
alently y) direction for 7 = 0.04 and 0.0025. The collapse is 
very good for the lowest shear rate. 
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for 7 = 0.04 and 0.0025. The collapse is very good for 0.0025 
only, while there is a systematic evolution for 0.04. 



VI. DISCUSSION 

We have investigated in this paper the ordering kinetics 
of a binary mixture quenched below its spinodal line in a 
homogeneous shear flow, through the numerical solution 
of the Cahn-Hilliard equation in 2 spatial dimensions. 

We have found that the bicontinuous coarsening struc- 
ture is well described by three parameters. Since the 
average shape of the domains is elliptic, it is sufficient 
to give the lengths Rj_{t) and i?||(<) of the two principal 
axis of the ellipse, and the orientation of the large axis 
with respect to the flow, 0{t). We have measured the 
time evolution of these three quantities, defined from the 
two-point correlation function, and all our results may be 
summarized by the following relations, which have been 
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shown to hold in the asymptotic regime of large times 

R\\{t)-it, 

R±{t) ~ constant, 



C(r,<) ~C 



(14) 



r|| rj_ 



In particular, this asymptotic regime could not be 
reached for the largest shear rates we have investigated. 
Moreover, our attempts at detecting (in the cases 7 = 
0.04 and 0.02) a crossover from a preasymptotic regime 
with a|| > 1 and a± > towards the regime (^4|) were 
unsuccessful, because this would require too large a sys- 
tem size Lx- These results also demonstrate that the 
domain growth does not stop under the shear flow, in 
the case where hydrodynamics is neglected and answer 
then questions (l)-(4) of the introduction. 

We now compare our results with previous ones. It is 
interesting to note that, even though this simulation was 
not intended to be able to reproduce experiments since 
hydrodynamics has been neglected, our results reproduce 
quantitatively the experiments of Chan et al. [||jl^ on a 
mixture of water and isobutyric acid and of Qiu et al. ||2l| ] 
on a polymer blend of polystyrene and poly- (vinyl methyl 
ether), for the three quantities describing the morphol- 
ogy of the domains. These experiments are probably 
done in a regime where hydrodynamic effects are neg- 
ligible, since other experiments by Hashimoto and col- 
laborators report a saturation of the domains to a 
7-dependent size. This stationary steady state has been 
termed 'string phase' and has also been observed by Hob- 
ble et al. 1^ in a strong shear rate regime. 

Our numerical findings can be compared with previ- 
ous works neglecting hydrodynamics, and our results for 
the growth laws are similar to those found by Qiu et 
al. and Corberi et al. ||l^. Note that in this last 
referenc e |l2] , the growth laws are compatible with the 
regime (|14|), although they are compared with the be- 
havior R^{t) - t^/^ and Ry{t) ~ t^/'^. The authors 
admit however that the latter regime was not reached 
within the numerical time-window. Earlier simulations 
were not able to be very quantitative [p^ , po|j2^ , |2^ , al- 
though Padilla and Toxvaerd have suggested that the 
anisotropic domain growth was well-described by alge- 
braic laws. This is, to our knowledge, the first time that 
the elliptic average shape of the domains is systematically 
investigated in a numerical experiment. Concerning the 
sca ling of the the two-point correlation function, Qiu et 
al. pifl have shown that the form (|l2| ) is not appropriate, 
but did not investigate the form (|13|). 

Analytically, the known results in the case of a con- 
served order parameter are for the 0{n) model, which 
has been solved with shear in the large-n limit |2J] (see 
also ||l^). The second result stems from a scaling argu- 
ment developed by Corberi et al. [n2| . Both predict that 



the growth laws should be Ra,{t) ~ t'^/'^ and Ry{t) ~ t^/^, 
with a scaling form for C{r,t) as in Eq. (|l2|). Our results 
do not corroborate these predictions. A reason may be 
that these results do not take into account the elliptic 
shape of the domains and we have clearly demonstrated 
this could be a key point for understanding the scaling 
properties of the system. One might also question the 
validity of an RG-type argument in the case where the 
exponent a± is zero, since upon space rescaling, domains 
become thinner and thinner. Recently psf , an approx- 
imate solution of the case of a non-conserved order pa- 
rameter has been found, predicting the scaling (13) in 
two spatial dimensions. However, extending this solu- 
tion to the conserved case is certainly a very hard task, 
since already in the unsheared case the conservation law 
makes the calculations very involved Q . 

We are then faced with the problem of having numeri- 
cal results that cannot be understood within an existing 
framework. We now give a simple physical argument 
leading to Eqs. (^4|) , inspired by the original argument 
given by Huse [pl|| to describe the zero-shear case. The 
spinodal decomposition under a shear flow results basi- 
cally from two competing effects. 

(i) The advection of the order parameter, which be- 
comes efficient for times t > 7^^, deforms the nearly 
circular domains existing at time t ~ 7"^. It is very 
easy to compute that, with the advection only, a circular 
domain is deformed into an elliptic shape with princi- 
pal axes scaling at large strain as R±{t) ~ (7^)^^ and 
i?l|(i) ~ 7t, with a tilt angle e{t) ~ (7^)"^ 

(ii) The domain growth arises because of the existence 
of a gradient of the chemical potential fj, = SF/6(j). This 
force gives rise to currents which make an interface of 
curvature R move with a velocity di?/dt cx 1/i?^. This 
interface motion results, in the absence of shear, to a 
coarsening of the domain structure Here, we mod- 
ify the argument by taking into account the fact that 
the structure is no longer isotropic. While the domain 
growth does not affect the tilt angle 0{t), it leads to two 
different interface velocities, namely di?||/di ^ — l/i?j_^ 
and dR±/dt ~ +l/i?||^ 

With the strong assumption that a balance can be 
made between these two effects, this leads to the follow- 
ing equations for the the three parameters of the ellipse 



dR^{t) 
dt 

dt 

e{t) 



1 



1 



7i2 i?||(i)2' 



= +7- 
1 



Rdt?' 



(15) 
(16) 
(17) 



which indeed imply Eqs. (|Tj). Although simple and 
heuristic, this argument captures in fact the essence of 
the coarsening process under a shear flow. It should cor- 
rectly describe the domain growth when the domains are 
quite large. Otherwise, the notion of interface velocity is 
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meaningless. It is quite clear however that subtleties, 
such as logarithmic corrections that exits in more in- 
volved computations , will not be captured by such a 
naive argument. 

Let us note that when it is applied to the non-conserved 
case, this argument leads to the growth laws R± ^ 
constant and i?|| t, at leading order. This not so far 
from the result R±_ ~ (Ini)-!/^ and t(lni)i/'' found 
in Ref. |2|. 

The second point we want to notice is that these equa- 
tions may be used to understand the following experi- 
ment. We have mentioned above the possibility of ap- 
plying the shear flow a certain time after the quench, in 
order to have lar ge c ircular domains as initial conditions. 
From Eqs. (|l5|)-(p7|), we expect the following behavior. 
For a relatively small strain, since the domains are quite 
large, the second term in the rhs of Eqs. (|l^) is negligible 
and Rj^ first decreases. Thus the second term in the rhs 
of Eq. becomes important and the growth of _R|| is 
slower than for a direct quench in the flow. This initial 
behavior is qualitatively different from the previous case 
and after this unusual transient, the growth becomes sim- 
ilar to the one studied previously. This predicted behav- 
ior agrees satisfactorily with the corresponding numerical 
experiment, as can be seen in Fig. Note, in particu- 
lar, that the slope of the curve i?|| (t) is minimum when 
R±{t) reaches its smallest value, as can also be deduced 
from Eq. (16). 
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FIG. 13. Time evolution of the two length scales as a func- 
tion of the strain when the system is first allowed to grow 
without shear during a time t = 400 and then submitted to 
a shear flow with 7 = 0.04. We have added as a visual guide 
a horizontal dashed line and the data for (t) for a direct 
quench in the shear flow (lowest curve). 



In conclusion, we hope this work has clarified several 
issues concerning the spinodal decomposition in a shear 
flow. It would be very interesting to perform the same 
measurements in three dimensions, since the analytical 
approach of Ref. psl shows that dynamic scaling proper- 



ties might be different in 2 and 3 dimensions. Of course, 
the next (big) problem is to have a better understanding 
of the effect of hydrodynamics on the phase separation 
in a shear flow. 
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